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Generating admissible space-time meshes for moving domains 
in d-\- 1-dimensions 

Elias Karabelas • Martin Neumuller 


Abstract In this paper we present a discontinuous Ga- 
lerkin finite element method for the solution of the 
transient Stokes equations on moving domains. For the 
discretization we use an interior penalty Galerkin ap¬ 
proach in space, and an upwind technique in time. The 
method is based on a decomposition of the space-time 
cylinder into finite elements. Our focus lies on three- 
dimensional moving geometries, thus we need to trian¬ 
gulate four dimensional objects. For this we will present 
an algorithm to generate d + 1-dimensional simplex 
space-time meshes and we show under natural assump¬ 
tions that the resulting space-time meshes are admis¬ 
sible. Further we will show how one can generate a four¬ 
dimensional object resolving the domain movement. First 
numerical results for the transient Stokes equations on 
triangulations generated with the newly developed mesh¬ 
ing algorithm are presented. 
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1 Introduction 

The finite element approximation of transient partial 
differential equations is in most cases based on explicit 
or implicit time discretization schemes. In particular 
the simultaneous consideration of different time steps 
requires an appropriate interpolation to couple the so¬ 
lutions at different time levels. Especially for spatial do¬ 
mains with a moving boundary one encounters various 
numerical difficulties. One usually relies on an arbitrary 
Lagrangian-Eulerian formulation. See for example the 
recent article m and references therein for an overview 
of the ongoing discussion. In this paper we consider the 
application of finite elements in the whole space-time 
cylinder Q. By this we mean a decomposition of Q into 
simplical elements. Therefore one replaces the problem 
of time discretization by a meshing problem. Having 
this, one can resolve the possible movement of the do¬ 
main Q directly. Simplicial space-time meshes have ad¬ 
vantages over tensor-product meshes, since it is easier 
to decompose complex space-time meshes by those ele¬ 
ments. 

Space-time finite element methods have been ap¬ 
plied to several parabolic model problems. Least square 
methods for convection-diffusion problems are consid¬ 
ered, e.g., in m and for flow problems, e.g., in m 
281122 j . Discontinuous Galerkin finite element 
methods have been applied to solve transient convection- 
diffusion problems in [24], for fluid dynamics see m 
and problems from solid mechanics see pii^irniiT] . Very 
recently a paper concerning the generation of 4D sim- 
plicial meshes from a sequence of 3D MRI data has 
been considered in [12]. Also rather recently, the X- 
FEM method has been considered in the space-time 
setting see mi. In most cases, the time dependent equa¬ 
tion is discretized in the space-time domain on space- 
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time slabs. This allows for local mesh refinement in the 
space-time domain, see for example EH- 

In this paper we consider, similar to [7], a decom¬ 
position of the space-time cylinder into simplicial finite 
elements. In particular for spatial domains i? C [R^ the 
space-time cylinder is a four-dimensional object, which 
has to be decomposed into finite elements. 

In [7], a method based on Delauny’s algorithm is 
given to construct a four-dimensional triangulation out 
of a given decomposition of the spatial domain i?. This 
method relies on the extension of the given finite el¬ 
ements of the triangulation of f2 to four-dimensional 
prisms. Afterwards a random perturbation of the re¬ 
sulting points is made, to ensure the admissibility of 
the resulting four dimensional mesh. Here, we present 
a different approach using similar ideas. Our method 
does not rely on random perturbations. Furthermore we 
can ensure and proof that the resulting mesh is admis¬ 
sible and we can also include movements of the domain 
boundary. We want to stress out, that our approach is 
still limited to a special class of boundary movements 
which we will describe in Section [3.41 

We will consider Stokes flow as a motivating model 
problem. For the approximation of the transient Stokes 
equations in the space-time cylinder we consider a dis¬ 
continuous Galerkin finite element method. In particu¬ 
lar, we apply an interior penalty approach in space [3l 
I6ll23[l21j . and an upwind technique in time 

This paper is organized as follows. In Section we 
describe the discontinuous Galerkin finite element method 
to solve the transient Stokes equations as a model prob¬ 
lem. The core part of this paper and the main results 
are given in Section where we describe our algorithm 
to generate a four-dimensional triangulation out of a 
given three-dimensional one. In Section we present 
some numerical results which underline the applicabil¬ 
ity of the proposed approach. We close the paper with 
some conclusions and comments on further work. 


2 Space-time discontinuous Galerkin Method 


Further we define the space-time mantle U as 
r ;= {{x,t) e ; a; e r{t) t e (0,T)} . 
The decomposition 0 induces 
U = Ud U Ur. 


The model problem we intend to study is governed by 
the transient Stokes equations. It reads as find (it,p) 
such that 


d 

—u - uAu ^\/p = f 
div(n) = 0 

u = gD 

Vn • n + aRU — pn = 

u = Uq 


in Q, 

in Q, 
on Ur, 
on Ur, 

on Uq := i?(0). 


( 2 ) 


Remark 1 In the case of a non-moving domain the def¬ 
inition of Q and U simplifies to 

Q:=nx (0,T), 

U :=dnx (0,T). 

For deriving a discrete variational formulation we need 
to decompose the space-time cylinder Q into simplicial 
elements, see [20]. Let 7^ be a sequence of decomposi¬ 
tions 

N 

Q = Th= IJ Tfc 

k=l 

into finite elements of mesh size hk- For d = 1 we have 
triangles, for d = 2 we use tetrahedrons and for d = 3 
pentatopes are chosen. The generation of such triangu¬ 
lations from a given triangulation of i7(0) is not trivial. 
We will address this topic in Section 

Definition 1 (Admissible decomposition) A de¬ 
composition Th is called admissible iff the non-empty 
intersection of the closure of two finite elements is ei¬ 
ther an edge (for d = 1,2,3), a triangle (for d = 2,3) 
or a tetrahedron (for d = 3). 


For any t G (0,T) let i?(t) C [R^ with d = 1,2,3 be 
a bounded Lipschitz domain with boundary r{t) := 
df2{t). We assume that the boundary r{t) admits the 
following decomposition for every t G (0, T) 

r(t) = r^(t)ur^(t). (i) 

We assume that the movement of the domain i?(t) is 
known for every t G [0,T]. We define the space-time 
cylinder Q as 

Q := G [R^^^ : x G i7(t) t G (0,T)} . 


It is worth noting that discontinuous Galerkin methods 
are not restricted to admissible decompositions. How¬ 
ever one needs additional technical assumptions, see |9]. 

Definition 2 (Interior facet) Let 77 be a decompo¬ 
sition of Q into finite elements For two neighboring 
elements e Th call 

Tki •= Tk n Ti 

an interior faeet iff r^i is a d-dimensional manifold. The 
set of all interior facets will be defined as Xh. 
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Any interior element has an exterior normal vector 
riki with a non-unique direction. We fix the direction of 
the normal vector such that riki is the exterior normal 
vector of the element Tk when k <1. So the direction of 
the normal vector riki depends on the ordering of the 
finite elements, but the variational formulation which 
we are going to use will be independent of this ordering. 


Definition 3 Let G be an interior facet with 
outer normal rik = G [Rd+l 

Tk and 

Til = —rik for Ti. For a given function (j) smooth enough 
restricted to either Tk or ti one defines : 

— The jump across Fki as 


The individual components read as 

N 


ah{uh,Vh) dq 

1=1 

/ {^x^h)Pkih^hirki^x 
. Jni 


and 




— V 


rki,x ^'^<75 




Fkieih 


+ 


/ 0:R{x,t)Uh 
JSr 


Vh dSq, 


l<Hkl ■= (t>\rknk + (l)\rini. 

— The space jump across Fki as 

•“ 0 |Tfc^x,/c + ^Irz^ccd* 

— The time jump across Fki as 

•“ 4^\rk'^t,k + 0|rz^t,Z' 

— The average of j) on Fki as 

i^)kl 2 ^\n'^ • 

— The upwind in time direction of (j) is defined as 




>|r. 

< 0 

<^\n 


if nk,t > 0 
if rik,t = 0 
if nk,t < 0 


Let p, g G Nq. Then one defines the spaces of piecewise 
polynomials 


■■= [SlirH)f = {vH e [LHQ)f : Vk\n € iWf 

for all n eTh, Vh\ED = 0 ] ^ 

Ql := G L^iQ) : qh\Ti ^ Pq{n) for all n eThj ■ 


bT{uh,Vh):='T- Uh-^Vhdq+ Uh-Vhdsg 

I—I dri J Et 

+ F [ dSq 

rkieih 

for a given velocity stabilization parameter au > 0. Fur¬ 
thermore we define the following pressure bilinear forms 
for Vh G V^{Th) and {ph,qh) G Ql{Th) x Ql{Th) : 


N 

bp{vh,Ph) •= T] / Phdiy{vh) dq 

;=i dri 

-F iPh)rJ^hlr,,,a.dsq, 

Fki 

dp{ph,qh)-= Yi Fki lPhirki,xMrk,,x dsq 

Fkieih drki 


for a given pressure stabilization parameter ap. In all 
the bilinear forms defined above we have used hki := 

^{hk + hi). Hence we have to find G V^{Th) and 

Ph G QliTh) such that 

A{uh,Vh) - bp{vh,Ph) = {f,Vh)Q + {uo,Vh) (3) 

-A{u'^,Vh), (4) 

bp{ul, qh) + dp{ph, qh) = -bp^, qh)- (5) 

Here we used an discrete extension of the given 
Dirichlet data, for example a i^-projection. 


3 Triangulations in d + 1 dimensions 


Inspired by works in min] we will use the following 
bilinear form defined for Uh^Vh G V^{Th)- 

A{uh, Vh) ■■= briuh, Vh) + ah{uh, Vu). 


In this section we will introduce an algorithm to de¬ 
compose a hyperprism into simplices to generate a (i+1 
simplex space-time mesh. Moreover we will show that 
the resulting mesh is admissible if the nodes of the sim¬ 
plices from the initial mesh are sorted in a special way. 
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3.1 Tensor product extensions 

A simple idea for constructing a space-time mesh for 
a given three-dimensional simplicial spatial mesh is to 
extrude the mesh in time direction by a tensor product 
extension, see also Figure Afterwards we decompose 
the upcoming prisms or so called hyperprisms into sim- 
plicial elements. 



Fig. 1: Tensor extension of a two-dimensional simplex 

Before we can start we need a precise definition of 
a d-dimensional simplex. 

Definition 4 (d-dimensional simplex) Let 
{Pi,---,Pd+i} c 

d G N be a set of nodes, then a d-dimensional simplex 
Sd is defined as 

Sd := [Pi, • • ■,Pd+i\ ■= conv({pi,.. .,Pd+i}), 

where conv(') is the convex hull of a set of nodes. Note 
that we also fix the ordering of the nodes in the defini¬ 
tion of a d-dimensional simplex. 

Now we can extrude one simplex in time direction 
and we obtain the following definition. 

Definition 5 (Hyperprism) For a given simplex Sd = 
[Pi^ ... ,Pd-\-i] the tensor product extension in time di¬ 
rection for a given time interval [0, r] or the so called 
hyperprism is given by 

Pd-\-l • [Pl •)'•'•) Pd-\-l 5 

:= conv({jj'i,... ... ,Pd+i}) C 

with 

p'i ■= {pI,oV, 

p" := 

for i = 1,..., d + 1. 


3.2 Decomposing Hyperprisms 

In this section we will give an algorithm to decompose 
the hyperprisms given in Definition into simplices. 

Definition 6 (Decomposed hyperprism) Let Sd be 

a given simplex and Pd-\-i the hyperprism with respect 
to the simplex Sd and r > 0. Then we define the fol¬ 
lowing simplices 

Sl+l ■= [pi,P2-P3>--->Pd+l-Pl]> 

S^+1 ■■= [P2,P3> • • • ,Pd+l,Pl,P2], 

•= [Ps; • • • ; P'd+1 ; Pl > P2 > Ps ]; (5) 

'S'dtJ := [Pd+ 1 , Pl, P2, P3, • • • , Pd+l]• 

Furthermore we define the set of simplices Tp{Sd, t) := 

{51+1,..., 

Note, that the ordering of the nodes of a hyperpism 
Pd+i is essential for the resulting decomposition 
In order to ensure that the simplices ..., 5'^^ 

defined in ^ decompose the hyperpism Pd-\-i we need 
the following lemma. 

Lemma 1 Let Pd-\-l be some given hyperprism with re- 
speet to the simplex Sd and r > 0. Then the set of sim- 
pliees 

rp{Sd,T) = {sl+„...,s^,Xl} 

defined in & is an admissible deeomposition of the hy¬ 
perprism Pdpi- 

Proof By construction the set of simplices Tp{Sd,T) = 

• • •, } is admissible. Furthermore, every sim¬ 

plex 5'^^! for i = 1,..., d + 1 is contained in the hyper¬ 
prism Pdpi since P^+i is convex. It remains to show, 
that the union of all simplices Tp{Sd^ t) is equal to the 
hyperprism, i.e. we have to show, that the volume of 
the union of all simplices Tp{Sd,T) coincides with the 
volume of the hyperprism. To do so, we transform the 
hyperprism Pdpi to a reference hyperprism Pdpi where 
we easily can compute all the volume terms. For this, 
we define the reference Simplex Sd C as 

Sd := [eo,ei,...,ed] = conv({eo, ei,..., 6^}), 

with 

eo := (0,0,...,0,0)T, 
ei ;=(1,0,...,0,0)T, 

62 := (0,1,...,0,0)T, 

6^:= (0,0,...,0,1)T. 
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Then we define the reference hyperprism Pd-\-i as 


Definition 7 (Consistently numbered) Let 


d+l 


[eo, • • •, !]• 


With the standard affine transformation we have a bi- 
jective mapping between the reference hyperprism Pd-\-i 
and the hyperprism This affine transformation 

consists of the standard transformation for d-dimensional 
simplices and a scaling in time direction. So we only 
have to compare the volume for the reference hyper¬ 
prism. Now the volume of the reference simplex Sd is 


given by 




= Hence the volume of the reference 


hyperprism is 


Pd+l 


1 


d\' 


The simplices of our decomposition in the reference do¬ 
main are given by 


Sl+i := 

[ cq , e'l, 62, . 

• • 5 ^0] 


[ei,e2,.... 

^d^ ^05 


[^25 • • • 5 

p" p" p" 
^05 5 ^2. 






It is easy to see, that these simplices have the same 
volume, i.e. 



1 

(dTi)!’ 


for i = 1 ,..., d + 1 . 


Hence we have 


d -\-1 




(d + 1 ) 


1 

(dTi)! 


1 


Pd+l 5 


which completes the proof. 


3.3 Admissible tensor product triangulations 

For a given d-dimensional triangulation Th we now want 
to construct a tensor product extension by applying the 
algorithm <© for every simplex of the simplicial mesh 
Th- With Lemma we know, that every hyperprism 
can be decomposed admissible into simplex elements. 
In this section we want to formulate conditions such 
that the overall space-time mesh is admissible. For this 
we need that the nodes of the simplices are ordered in 
a special way. 


be an admissible d-dimensional simplex mesh. Then Th 
is called consistently numbered^ iff for any two simplices 
^ Th with non-empty intersection, i.e. >5^ fl >5^ 7 ^ 
0, there exists indicies ki < ... < kn and < ... < 
with n G N, n < d + 1 , such that 

5'd n .. .,pIJ = \p\^,... ,pIJ. 

Here means that the two sets are the same and 

means that the two sets are equal and that also 
the numbering of the nodes is the same, i.e. = 

The definition of a consistently numbered triangulation 
can also be found in [ 8 ] and it is important for the 
refinement of d-dimensional simplices, especially for d > 
4. If an admissible mesh is consistently numbered we 
can prove the next Theorem. 

Theorem 1 Let Th be an admissible d-dimensional tri¬ 
angulation which is consistently numbered and let r > 
0. Furthermore let 

Th^r-={Tp{Sd^T):SdeTh} 

be the (d + 1)-dimensional simplex mesh resulting by 
decomposing every hyperprism with the algorithm given 
in a. Then the space-time mesh Th,r 'Is admissible. 

Proof With Lemma [2 we know, that every hyperprism 
is decomposed admissible into simplices. To obtain a 
global admissible mesh we have to prove, that the tensor 
product triangulations Tp(5'^, r) and Tp(S'^, r) for each 
neighboring elements Sp^d^Ph are matching. Let 

Sd = \Pi,---,Pd+i] and 5^ = 

with 5'^, 5'^ G 7 ^ be some neighboring simplices and 

Pd ■■= Prf+i 

with 

Pd+i ■= [Pi> • • • >Pd+i; 't] and Pj_^i := \p{,.. • ,p^+i; t] 
be the intersecting hyperprism and 

n := {5d+i n Pd : Sd+I e TpiS^T)}, 
r/ := {Sd+i n Pd : Sd+i € Tp{Si,T)}, 

be the corresponding triangulations of Pd obtained by 
rp(5*,r) and Tp{S^^ , r). It remains to show, that the 
intersecting hyperprism Pd is decomposed in the same 
way from both sides, i.e. that Tp = Tp- Since Th is 
















6 


Elias Karabelas, Martin Neumiiller 


consistently numbered, there exists indices ki < ... < 
kn and < ... < In with n = such that 

~ ^d-l [Pki^ • • • ^PkJ 

='S'Li ;= (^) 

Therefore, the intersecting simplex fl 5'^ is obtained 
by simply removing the nodes from >5^ or which 
are not shared together and furthermore they have the 
same ordering of the nodes. For the intersecting hyper- 
pism Pd the decompositions from both sides Tp and Tp 
are given by removing the nodes which are not shared 
together from the formula (§ and with 0 we have 

n = Tp{SU,r) and = Tp{Si_„T). 

Since in equation Q also the node ordering of 
and is the same we also obtain 

which implies that Tp = Tp. 

Remark 2 To obtain an admissible space-time mesh Th^r 
we only have to ensure, that the nodes of the spatial 
mesh Th are consistently numbered. This can be easily 
obtained by sorting for each simplex Sd ^ Th the local 
nodes with respect to the global node numbers. 

3.4 Tensor product triangulations for moving domains 

If the movement of a computational domain is known in 
advance we can generate admissible space-time meshes 
by applying the methods from above. The idea is to 
move the points at the top of the tensor-product exten¬ 
sion. Assuming that the displacement of points on the 
boundary T(t) is governed by a function 

g^,dX,t):r{0)x{0,T)^R^. 

Then a point x e r{t) can be written as 

X X T ^ t), 

where X G T(0). Recall the definition of a hyper¬ 
prism in Definition]^ Instead of using p'l := (p7,r)^ 
on the surface we can apply the displacement and use 
p'l := {pjPOmoY^Pi^ all boundary points of 

the simplex mesh that are subject to a movement. The 
remaining generation of the 4D mesh stays untouched. 
For boundary movements that are of small magnitude 
and do not change the topology of the initial geometry 
this can be sufficient. For stronger yet topology preserv¬ 
ing movements this concept would create degenerating 
simplex elements. A remedy to this is to use the move¬ 
ment S^rnov as Dirichlet datum for a vector Laplacian or 


a linear elasticity problem. Then the resulting displace¬ 
ment is applied to all simplex points in the domain. For 
more on mesh smoothing we refer to mils]. 

In the case of stronger displacements or even topol¬ 
ogy changes re-meshing would be required and we need 
further meshing algorithms to connect different spatial 
domains in space and time, especially for four-dimensional 
space-time meshes this remains a future research topic. 

If the the movement of the computational domain 
is not known in advance we can solve the problem on 
a coarse spatial grid with coarse time steps to obtain 
a coarse approximation for the movement. Afterwards 
we can construct the coarse space-time mesh with the 
methods given in this work. By using adaptive schemes 
in space and time we further can refine the space-time 
domain adaptively and move the points in the space 
time domain by the computed finer approximations. 
Note that the movement of the points has to be only 
done in the range of the approximation error, which 
is usually small. Of course this is also considered as a 
further research topic. 

3.5 Visualization 

Here we want to address the issue of visualizing results 
for four-dimensional triangulations 7^. In applications 
it is desired to visualize results at given time instances 
tk G [0,T]. The main idea is to cut the decomposition 
Th into a finite number of three-dimensional manifolds. 
For this we need to have a hyperplane to calculate the 
intersections with the decomposition. 

Definition 8 (Hyperplane) Let Pq G be arbitrary 
and let P 2 , p^ and p^ G be linear independent. 
Then the set 

:= {x e R‘^ : X = Pq + + 112P2 + M3P3 

for ij.i,P2,P3 e R} 

is called a hyperplane. 

To cut a given decomposition Th with a hyperplane iL 4 , 
we have to cut every element Tk ^ Th with the hy¬ 
perplane. For this we have to calculate for every edge 
ei = {xi-^^Xi^)^ i = 1,..., 10 of T/c, the intersection with 
the hyperplane. A point x e e^ can be written as 

X — ^ii T A (^22 

for a given A G [0,1]. Hence, an intersection point of 
the edge e^ with the hyperplane H 4 has to satisfy 

Xi, + A {Xi^ -Xid =Pq + PiPi + P2P2 + M3P3 
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or in matrix notation 


4.1 Robin Boundary Conditions for Simulating Valves 


Aj ;= (pi P2 Ps Xi 


Xj. 


(Pl\ 

P 2 

Va/ 


Xj 


■Po- 


The matrix is invertible iff the vector Xj^ — Xj^ is 
linear independent to the vectors P 2 , P 3 . In fact, 
the matrix is not invertible if the edge Cj is parallel 
to the hyperplane In this case there exists either 
no intersection point or infinitely many. If the matrix is 
invertible we can calculate the coefficients and 

A G [R uniquely. Let denote the set of all intersection 
points of the element ^ Th with the hyperplane 
We distinguish two relevant cases 

1. If \Dk\ = 4, then the intersection points form a 
tetrahedron 

2 . If I = 6 , then the intersection points for a general 
irregular prism. 

If we use the special vectors 


Po •— Pi •— ^x-)P2 •— 

for a given t* G [0, T] we can now calculate a three- 
dimensional object which can be visualized with exist¬ 
ing software tools for example [4]. 


In the subsequent examples we want to simulate open¬ 
ing and closing valves. This means, that we need to 
switch between and inflow and a no-slip condition. To 
this end we used the following configuration for §: We 
set 9r = 0- Further we use the following Robin coeffi¬ 
cients for outflow 


aR{x,t) 


10 ® ifiG[0, i) 

0 if t e 1 ] 


and the following for inflow: 


aR{x,t) 


0 if f G [0, |) 
10® iffG[i,l]' 


4.2 First Example 

In the first example we consider Stokes flow in a di¬ 
aphragm pump. The geometry consists of the intersec¬ 
tion of two cylinders. The first one has its main axis 
aligned with the z-diX.is with a radius of 0.8 and ranges 
between 2 : = —0.4 and z = 0.4. The second cylinder has 
its main axis aligned with the x-axis with a radius of 
0.2 and ranges from x = —ltox = l.A front view of 
the geometry is depicted in Figure The movement of 


4 Numerical Results 


In this section we will present first numerical examples. 
Starting point is the discrete variational formulation 
(i-©. This can be equivalently written as the follow¬ 
ing block system 


VB/^ J 



(8) 


It is worth noting, that due to the discretization of the 
time derivative we have that 7 ^ K^. The four dimen¬ 
sional computational geometries as well as the resulting 
linear systems in the subsequent numerical examples 
were solved with the software package Neshmet devel¬ 
oped by the authors. In particular we used a precondi¬ 
tioned GMRes method. As preconditioner we used the 
following: 




Fr 



rR,i 


Fig. 2 : Front view of the initial geometry i?(0). Blue 
boundaries belong to Fd- 


the boundary r{t) was prescribed as follows 


9mov{t,X) := 0.4 + sin^TTf) 1 - 


V(0)2 + X(l)2 


0.752 


-X 


P := 



for X G TD,m(0) and 0 else. The following boundary 
conditions are used: 


where is chosen as a component-wise algebraic multi¬ 
grid and Sh is chosen as ILU(2)-factorization of + 
B/j,diag(K/j,)“^B^. These preconditioners were taken from 
the HYPRE library m- 


u = 0 on Fd 
U = WtOmov on rD,m{t) 

On and TR,out we used the Robin boundary 

conditions discussed in Section ItT] 
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— The initial condition for u was set to n(0, x) = 0. 

The triangulation of the resulting 4D geometry was ac¬ 
complished using the tools described in Section The 
resulting mesh consisted of 951360 pentatopes. 


needed 95 GMRes-iterations for achieving a relative er¬ 
ror of — 5. In Figureone can see the resulting flow 
and pressure at given time stamps, which were again 
produced by slicing the 4D geometry along the time 
axis. 



(c) (d) 

Fig. 3: Initial Triangulation i?(0) 


Some snapshots of the triangulation of the moving 
domain are depicted in Figure]^ This snapshots were 
generated by slicing through the 4D mesh along the 
time axis as described in Section [Tb) The polynomial 


4.3 Second Example 


For the second example we considered a Y-shaped pipe. 
A schematic view is depicted in FigureWe prescribed 
the following movement of T(t): 


9m.ov^^I • 


l^(2)+3| 


sin(7rt)^e^ 


for X ^ U 
for X G Tc),m U Fb^t, 


Some snapshots of the domain movement are depicted 
in Figure The boundary conditions were set to 

— n = 0 on F]j U FD,m{t) 

— u = on 

— On and TR,out we used the Robin boundary 

conditions discussed in Section ItT] 

— The initial condition for u was set to u{0^x) = 0. 

The resulting 4D mesh consisted of 2618880 pentatopes. 
With the same ansatz spaces as used for Example 1 we 
have 39283200 degrees of freedom for and 2618880 
degrees of freedom for We needed 107 GMRes-iterations 
for achieving a relative error of IE — 5. In Eigurej^we 
have depicted some results. 



(a) 12(0.5) (b) 12(0.5) 



(c) 12(0.65) (d) 12(0.8) 

Eig. 4: Snapshots of the Triangulations 12(t) 

degree for Uh was set to p = 1 and q = Ofoi the pressure 
variable ph- This resulted in 14270400 degrees of free¬ 
dom for Uh and 951360 degrees of freedom for ph- We 


5 Conclusions 

In this paper we have presented a novel approach to 
construct four-dimensional triangulations for moving 
domains. This was done by extending the elements of 
the space triangulation into hyperprisms. Assuming a 
consistent numbered spatial triangulation we were able 
to proof that our algorithm produces admissible space- 
time meshes. We implemented the presented algorithm 
and applied it to solve the transient Stokes equations 
with a space-time discontinuous Galerkin finite element 
method. In the future one could start investigating Navier- 
Stokes equations. Eurthermore, optimal control prob¬ 
lems with time dependent partial differential equations 
render themselves interesting candidates for applying 
space-time methods, since one has to solve a forward 
and backward problem which are coupled in space and 
time. Another attractive aspect of general space-time 
meshes is the possibility to apply adaptive refinement 
strategies to resolve local behaviors in space and time. 
Gonsidering solvers, one could think about domain de¬ 
composition approaches or space-time multigrid meth¬ 
ods for example, which are a future research topic. 
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(c) t = 0.7 (d) t = 1.0 

Fig. 5: Snapshots of the solution. Additionally we have cut along the ^-axis. 
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Fig. 7: Initial Triangulation i7(0) 
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(a) 12(0.19) (b) 12(0.54) 

r Y 

(c) 12(0.77) (d) 12(0.95) 

Fig. 8; Snapshots of the Triangulations Q{t) 
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